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Abstract 

We study characteristic (quasinormal) modes of a .D-dimensional Schwarzshild 
black hole. It proves out that the real parts of the complex quasinormal modes, 
representing the real oscillation frequencies, are proportional to the product of 
the number of dimensions and inverse horizon radius ~ D r^ 1 . The asymptotic 
formula for large multipole number / and arbitrary D is derived. In addition the 
WKB formula for computing QN modes, developed to the 3th order beyond the 
eikonal approximation, is extended to the 6th order here. This gives us an accurate 
and economic way to compute quasinormal frequencies. 



1 Introduction 



Within the framework of the brane world models the size of extra spatial dimensions 
may be much larger than the Plank's length, and the fundamental quantum gravity scale 
may be very low (~ Tev). When considering models with large extra dimensions the 
black hole mass may be of order Tev., i.e. much smaller than the Plank's mass. There 
is a possibility of production of such mini black holes in particle collisions in colliders 
and in cosmic ray experiments pQ. Estimations show that these higher dimensional 
black holes can be described by classical solutions of vacuum Einstein equations. Thus 
the investigation of general properties of these black holes, including perturbations and 
decay of different fields around them, attracts considerable interest now (see for example 
[2], [2] and references therein). 

It is well-known that when perturbing black hole it undergoes damping oscillations 
which are characterized by some complex eigenvalues of the wave equations called quasi- 
normal frequencies. Their real parts represent the oscillation frequencies, while the imag- 
inary ones determine the damping rates of the modes. The quasinormal modes (QN) of 
black holes (BH's) depend only on a black hole parameters and not on a way in which 
they were excited. QN's are called, therefore, "footprints" of a black hole. Being a useful 
characteristic of black hole's dynamics, quasinormal modes are studied also within differ- 
ent contexts now: in Anti-de-Sitter/Conformal Field Theory (AdS/CFT) correspondence 
(see for example jlj- [TH] and references therein), because of the possibility of observing 
quasinormal ringing of astrophysical BH's (see ^H] for a review), when considering ther- 
modynamic properties of black holes in loop quantum gravity J7|-[2n], m the context of 
possible connection with critical collapse jl], 0, 12] [22]- 

Thus it would be interesting to know, from different grounds, what happen with QN 
spectrum a black hole living in D -dimensional space-time [22] > The present paper is 
two-fold: First we extend the WKB method of Schutz, Will and Iyer for computing QN 
modes from the 3th to the 6th order beyond the eikonal approximation (see Sec. II and 
Appendix I). In a lot of physical situations this allows us to compute the QNMs accurately 
and quickly without resorting to complicated numerical methods. In Appendix II QN 
modes of D = 4 Schwarzshild black hole induced by perturbations of different spin are 
obtained by the 6th order WKB formula, and compared with the numerical values and 3th 
order WKB values. Second, motivated by the above reasons, we apply the obtained WKB 
formula for finding of the scalar quasinormal modes of multi-dimensional Schwarzshild 
black hole (Sec. III). It proves out that the real parts of the quasinormal frequencies are 
proportional to the product Dr^ 1 , where r is the horizon radius, and D is the dimension 
of space-time. 

2 Sixth order WKB analysis 

First semi-analytical method for calculations of BH QNMs was apparently proposed by 
Bahram Mashhoon who used the Poschl- Teller potential to estimate the QN frequencies 
[21]. In [23] there was proposed a semi-analytical method for computing QNM's based on 
the WKB treatment. Then in [21] the first WKB order formula was extended to the third 
order beyond the eikonal approximation, and, afterwards, was frequently used in a lot 
of works (see for example [H] [2E], [2H]- [SO], [21]; [22]; [22]; [21] and references therein). 
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The accuracy of the 3th order WKB formula (see eq. (1.5) in |26j) is the better, the 
more multipole number / and the less overtone n. For the Schwarzshild BH the results 
practically coincide with accurate numerical results of Leaver [33] at I > 4 when being 
restricted by lower overtones for which I > n. For fewer multipoles, however, accuracy 
is worse, and may reaches 10 per cents at I = 0, n = 0. Numerical approach on 
contrary, is very accurate, but, dealing with numerical integration or systems of recurrence 
relations, is very cumbersome, and, often, require modification to be applied to different 
effective potentials. At the same time WKB approach lets us to obtain QNM's for a full 
range of parameters giving thereby some fields of work for intuition as to physical behavior 
of a system. Even though WKB formula gives the best accuracy at I > n, it includes 
the case of astrophysical black hole radiation where only lower overtones are significantly 
excited [3^]. Both advantages and deficiencies of the WKB approach motivated us to 
extend the existent 3th order WKB formula up to the 6th order. 

The perturbation equations of a black hole can be reduced to the Schrodinger wave- 
like equation: 

^ + Q{x)i,{x) = 0, (1) 

where "the potential" —Q(x) is constant at the event horizon (x = — oo) and at the 
infinity (x = +oo) and it rises to maximum at some intermediate x = xq. Consider 
radiation of a given frequency u incident on the black hole from infinity and let R(u) 
and T(uj) be the reflection and transmission amplitudes respectively. Extend R{uS) to 
the complex frequency plane such that Re(z) ^ 0, and T(z)/R(z) is regular. Then, the 
quasinormal modes correspond to the singularities of R(z). We have a direct analogy with 
the problem of scattering near the pick of the potential barrier in quantum mechanics, 
where uj 2 plays a role of energy, and the two turning points divide the space into three 
regions at which boundaries the corresponding solutions should be matched. 

To extend the 3th order WKB formula of we used the technique of Iyer and Will. 
We shall omit here the technicalities of this approach which are described in The 
only thing we should stress is that since the coefficients My, that connect amplitudes near 
the horizon with those at infinity, depend only on v (related to the overtone number n) 
they may be found to higher orders, simply by solving the interior (between the turning 
points) problem to higher orders. Thus there is no need to perform an explicit match 
of the solutions to WKB solutions in the exterior (outside turning points) regions to the 
same order. The result has the form: 

'2L - A 2 - A 3 - A 4 - A 5 - A 6 = n + \ , (2) 

where the correction terms A 4 , A 5 , A 6 can be found in the Appendix I. Note that A 4 
coincides with preliminary formula (A3) of i n proper designations. 

An alternative, pure algebraic approach to finding higher order WKB corrections was 
proposed by O.Zaslavskii jH7j, using a quantum anharmonic oscillator problem where 
WKB correction terms come from perturbation theory corrections to the potential an- 
harmonicity. 

Thus we have obtained an economic and accurate formula for straightforward calcu- 
lation of QNM frequencies. The 6th order formula applied to the D = 4 Schwarzshild 
BH is as accurate already at I = 1 as the 3th order formula does at / = 4. We show it in 
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Appendix II on example of QNM's corresponding to perturbations of fields of different 
spins: scalar (s = 0), neutrino (s = 1/2), electromagnetic (s = 1), gravitino (s = 3/2), 
and gravitational (s = 2). In addition, looking at the convergence of all sixth WKB 
values to some unknown true QN mode, we can judge, approximately, how far from the 
true QN value we are, staying within the framework of WKB method. 

3 Quasinormal modes of the D-dimensional 
Schwarzshild black hole 

The metric of the Schwarzshild black hole in D-dimensions has the form: 

ds 2 = f{f)dt 2 - f-\r)dr 2 + r 2 dQ 2 D _ 2 , (3) 

where 

f( \ 1 f r °^~ 3 1 167tGM U\ 

Here we used the quantities 

(2n)( D - 1 V 2 

nD - 2= T((D-l)/2Y r ( 1 / 2 ) = ^- r(z + l) = zT(z) 

The scalar perturbation equation of this black hole can be reduced to the Schrodinger 
wave-like equation (JTJ) with respect to the "tortoise" coordinate x: dx = jfr where "the 
potential" —Q(x) has the form: 



Q(x) =u 2 - f{r) ( ^ '- + ^ ^ '-f{r) + —^—f{ T ) ) , (5) 



At some fixed D we can put r = 2 and measure u in units 2r . The quasinormal modes 
satisfy the boundary conditions: 

<j>{x) rsj c±e ±iuJX as x-+ ±oo. (6) 

The 6th WKB order formula used here gives very accurate results for low overtones. 
The previous orders serve us to see the convergence of the WKB values of uo 2 as a WKB or- 
der grows to an accurate numerical result. Namely we can observe that for I = 1, 2, 3, 4, .. 
for the fundamental overtone the 6th order values differs form its 5th order value by frac- 
tions of a percent or less at not very large D (we are restricted here by D = 4, 5, ...15). 

It proves out that if one takes ro = 2 for each given D, then the real parts of uj 
for different D lay on a strict line. That is, UR e is proportional to the product r$D 
(Remember that r depends on D itself). Namely, for the fundamental overtone we 
obtain the following approximate relations: 

u Re ~ 0.244 J D(r /2)" 1 , 1 = 2 (7) 
u Re ~ 0.275D(r /2)-\ 1 = 3 (8) 
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uju^Q.290D{r /2)-\ 1 = 4. (9) 

Here we take uj = uj Re — iu>j m . Generally, the more the multipole number I, the more the 
coefficient before the product Dtq 1 . The same ~ Dr^ 1 relation we observed for higher 
overtone but not higher than I, for which WKB treatment is applicable. In Fig. 1,2 we 
presented the real and imaginary parts of u measured in 2r ( ^ 1 for different D. For real 
parts of I = 1 modes we see the deviation from the strict line at large D. This however, 
is stipulated by a bad accuracy of the WKB approach, and we believe that the true 
frequencies will lay on strict line again. Indeed, one can judge about it by looking at 
the convergence plot Fig.3-Fig.6 where the real and imaginary parts of u are shown as a 
function of the WKB order. Generally the accuracy of the WKB formula is the better, 
the more I, and the less n and D. Note that the dependence Dtq 1 for lower overtones 
can be recovered even within 3th order formula, provided I is greater than 2, and D is 
not very large. 

Another point is the I = modes: in this case the lowest overtone implies I = n, 
and the WKB formula has considerable relative error. For a four-dimensional BH, for 
which the accurate numerical results are known, the error is about 10 percent for u>j m , 
and 5 percent for uo Re in the third WKB order, while in the sixth order it reduces to 
percent for u Re and 3 percent for ujj m (see Appendix II). For greater D the error 
increases, the difference between the fifth and sixth order WKB values grows and one 
cannot judge of true quasinormal behavior in this case (see Fig.3, 4, 5.). Fortunately, 
other field perturbations, including gravitational, have the lowest overtone with I > n 
and the WKB treatment is of good accuracy for all I. In Table 1. we compare the third 
order WKB values of I = 0, n = modes for different D [3] with those obtained through 
the sixth order here. 
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Figure 3: uor c (bottom) and Ui m (top) as a function of WKB order of the formula with 
which it was obtained for I — 1, n — 2, D — A modes, and the corresponding numerical 
value. We see how the WKB values converge to an accurate numerical value as the WKB 
order increases. 
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WKBorder 



Figure 4: oo Re (top) and io Im (bottom) as a function of WKB order of the formula with 
which it was obtained for I = 0, n = 0, D = 12 modes. 




Figure 5: ujr c (top) and ui m (bottom) as a function of WKB order of the formula with 
which it was obtained for/ = 0,n = 0,-D = 6 modes. 



D 


3th WKB order 


6th WKB order 


l/r 


4 


0.1046 


- 0.11522 


0.1105 


-0.1008? 


0.5 


6 


1.0338 


- 0.7133i 


1.1808 


- 0.6438i 


1.28 


8 


1.9745 


- 1.0258i 


2.3004 


- 1.0328* 


1.32 


10 


2.7828 


- 1.1596i 


3.2214 


- 1.37662 


1.25 


12 


3.4892 


- 1.2020i 


3.9384 


- 1.75742 


1.17 



Table I. Schwarzshild QN frequencies for / = 0, n = scalar perturbations in various D. 
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Figure 6: uiR e (top) and uji m (bottom) as a function of WKB order of the formula with 
which it was obtained for I — 1, n = 0, D — 6 modes. 



For large I the well-known approximate formula reads (see j3H], [H] [HH] for a proof) 

"* = 3^s(' + 5) > 1 *" = 37S (" + D (10) 

To obtain its .D- dimensional generalization we find a value r max at which the effective 
potential V attains its maximum, provided I is large 



„ -D-4 

max ~ 



2^(5-1)5^, £> = 4,5,6,.... (11) 



Then let us make use of this value r max when dealing with the first order WKB formula. 
After expansion in terms of small values of 1/7, for a fixed D in units of 2tq 1 we obtain 



D + 2Z-3 / 2 \— /D-3 

(12) 



4 V-D-17 VD-1 



(D-3)/ 2 2n+l 

When D = 4 these formulas go over into pijl . We see that when I is much larger 
than D, the ~ Dr^ 1 dependence of tUR e breaks down. 



4 Conclusion 

We were interested here in a question how dimensionality effects on quasinormal behavior 
of black holes. Yet, several interesting points are left beyond our consideration of low 
laying quasinormal modes of multi-dimensional black holes. First of all, one would like 
to understand the origin of the relation ~ Dr^ 1 in ujR e dependence. In this question it 
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is possible to try to explain it from the interpretations of QN modes as Breit-Wigner 
type resonances generated by a family of surfaces wave propagating close to the unstable 
circular photon orbit jUJ]. Second, we do not know whether ~ Dr^ 1 dependence will be 
present for perturbations of other fields, and for more general backgrounds, such as multi- 
dimensional Reissner- Nordstrom or Kerr. We hope further investigations will clarify these 
points. 
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5 Appendix 1: Correction terms for WKB formula 

Here we shall follow the designations: Qo means the value of the potential Q at its pick, 
while Qi is the ith derivative of Q with respect to the tortoise coordinate x. Then Q\ is 
the jth power of the ith derivative of Q. 

A 4 = 1 _ - - (2536975Qg - 9886275Q 2 <9&4 + 5319720Q^Q?Q 5 - 

597196800^2^ 



225Q'g'(-40261Q^ + 9688Q 2 Q 6 ) + 3240Q 2 i Q 3 (-1889Q 4 Q5 + 220Q 2 Q 7 
729Q 3 2 {U25Ql - 1400Q 2 Q 4 Q 6 + 8Q 2 {-123Q 2 5 + 25Q 2 Q 8 )))+ 
(n + 1/2) 2 

: (348425g^ - 1199925g 2 g^g 4 + 57276QlQlQ 5 - 



4976640 v / 2g^^ 
A5QtQl{-20671Ql + A552Q 2 Q 6 ) + 1980Q 3 2 Q 3 {-A89Q A Q 5 + 52g 2 g 7 )- 

27g^(2845g^ - 236og 2 g 4 g 6 + 56g 2 (-3ig2 + sg 2 g 8 )))+ 

(n + 1/2) 4 

: (192925g^ - 581625g 2 g^g 4 + 234360g^g 5 - 



2488320 v / 2g^\/g^ 

45g2g2(-83i5g^ + i448g 2 g 6 ) + io8og3g 3 (-i6ig 4 g 5 + i2g 2 g 7 )- 

27g^(625g^ - 44og 2 g 4 g 6 + 8g 2 (-63g' + sg 2 g 8 ))) (14) 

(n + 1/2) 

A 5 = 573 3 089 2800gio (2768256Q 10 g 2 7 - 1078694575g« + 5357454900g 2 g^g 4 - 

2768587920g^g^g 5 + 90g^g^(-88333625g^ + 12760664g 2 g 6 )- 
4320g^g^(-1451425g 4 g 5 + 91928g 2 g 7 ) - 27g^(7628525g^ - 9382480g 2 g^g 6 + 
64g2(19277g2 + 37764g 5 g 7 ) + 576g 2 g 4 (-21577g2 + 2505g 2 g 8 )) + 
540g^g^(6515475g^ - 3324792g 2 g 4 g 6 + 16g 2 (-126468g^ + 12679g 2 g 8 ))- 

432g^g 3 (5597075g2g 5 - 854i6og 2 g 4 g 7 + 8g 2 (-i454i7g 5 g 6 + 6685g 2 g 9 )))+ 
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J^±^^(31lOAQ 10 Ql - 42944825Q« + 193106700Q 2 QiQ 4 - 
-90039120Q^Q 5 + 30g^Q^(-8476205g^ + 1102568g 2 g 6 )- 

432og3g3(-4ii65g 4 g 5 + 23i2g 2 g 7 ) - 9g 2 t (445825g* - 47288og 2 g;jg 6 + 
64g2(829g2 + i836g 5 g 7 ) + 4032g 2 g 4 (-i79gj? + isg 2 g 8 ))+ 
i8og^g2(5326i5g3 - 24i224g 2 g 4 g 6 + i6g 2 (-9352g2 + 799g 2 g 8 ))- 
i44g^g 3 (392325g2g 5 - 5i6oog 2 g 4 g 7 + 8g 2 (-8853g 5 g 6 + 335g 2 g 9 )))+ 

( n _|_ I/O) 5 

ii943936oog 2 o (10368gi ° g ^ " 66578225 ^3 + 272124300g 2 g 3 , g 4 - 

112336560g^g^g 5 + 9450g^(-33775g^ + 3656g 2 g 6 )- 

i5i2oog3g3(-i297g 4 g 5 + 56g 2 g 7 ) - 27g 2 1 (89075g* - 8344og 2 g;jg 6 + 

64g^(i3ig2 + 396g 5 g 7 ) + 576g 2 g 4 (-343g2 + isg 2 g 8 ))+ 
54og^g2(i88i25g^ - 7i4oog 2 g 4 g 6 + i6g 2 (-3052g2 + i77g 2 g 8 ))- 
432g^g 3 (ii8825g^g 5 - ii76og 2 g 4 g 7 + 8g 2 (-2303g 5 g 6 + 55g 2 g 9 ))) (15) 

Ae = 20226338979"8400g^ ( " 171460800gi2 ^ + 1714608000 ^^- 



10268596800g 10 g 2 g3 + 970010662775g^° + 3772137600g 10 g 2 g 4 - 
6262634175525g 2 g®g 4 + 13782983196150g^g^ - 11954148125850Q^Q|+ 
3U9l70577475Q 4 2 Q 2 3 Qi - 144528059025g^g^ + 3352602187200g^g^g 5 - 
12300730092000g^g^g 4 g 5 + 11994129604800g^g^g^g 5 - 2624788605600g 2 , g 3 glg 5 + 
2580769643760g^g^g^ - 3453909784416Q^Q 4 <95 + 438440697072g^g^g2 + 
+260524397952g 2 i g 3 gj? - 1475306441280g^g^g 6 + 4329682610400g 2 1 g 3 1 g 4 g 6 - 
2865128172480g^g^g2g 6 + 233443879200g^g 6 - 1660199804928g 2 , g3g 5 g 6 + 
1281705296256g 2 , g 3 g 4 g 5 g 6 - 87403857408g 2 7 gj?g 6 + 231105873600g^g2- 
68412859200g^g 4 g^ + 552968700480g^g 7 - 1231789749120g^g 4 g 7 + 
470726303040g^g 3 g^g 7 + 413953400448g^g^g 5 g 7 - 126242178048g^g 4 g 5 g 7 - 
91489305600g^g 3 g 6 g 7 + 5619715200g^g^ - 175752294480g^g 8 + 
271759652640g^g^g 4 g 8 - 39736040400g^g^g 8 - 73378363968g^g 3 g 5 g 8 + 
9773265600g^g 6 g 8 +47107126080g^g 9 -43345290240g^g 3 g 4 g 9 +7400248128g^g 5 g 9 )- 

(n + l/2) 2 i 

=(-4551552Qi 2 Q^ + 60279552QiiQ|Q 3 - 



687970713600g^ 2 v /2g 2 ~ 
425036160gi g 2 g 3 + 73727194625g^° + 116743680gi g 2 g 4 - 
443649208275g 2 g^g 4 + 901144103850Q&jjQ2 - 711096726150g^gl+ 
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182164306725Q4Q2Q4 _ 6289615575g 2 , g 4 ' + 2224Q7Q24W0Q 2 2 QlQ 5 - 
7AQA18AA5200QlQlQ 4 Q 5 + 653423900400g^g^g^g 5 - 124319674800g^g 3 glg 5 + 
143980943040g^g^ - 169712521920g 2 , g 3 ; g 4 g2 + 18188188416g^g2+ 
11240861184g^g 3 g3 - 91198200240g^g^g 6 + 241513732080g^g 4 g 6 - 
140030897040g^g2g 6 + 9200103120g^g 6 - 84218693760g^g^g 5 g 6 + 
55248386688g^g 3 g 4 g5g6 - 3173043456g^g2g 6 + 10464952896g^g^g2- 
2403421632g^g 4 g^ + 31637744640g^g^g 7 - 62649953280g^g3g 4 g 7 + 
20409822720g^g 3 g2g 7 + 18860532480g^g^g 5 g 7 - 4693344768g 2 7 g 4 g 5 g 7 - 
3625731072g^g 3 g 6 g 7 + 188054784g 2 i g2 - 9155635200g^g 8 + 

i223802432og^g2g 4 g 8 - i40527872og 2 / g;jg 8 - 28667ooi6og^g 3 g 5 g 8 + 

303295104g^g 6 g 8 + 2210705280g 2 i g 3 ; g 9 - 1685525760g 2 / g 3 g 4 g 9 + 235488384g^g 5 g 9 ) - 
( " + ' /2) ' ' -QQ528Q 12 Q 9 2 + 1245888QnQiQ 3 - 



200658 l2mQ¥^2Q~ 2 
11158560gi g2g 3 + 4668804525g^° + 2116800gi g2g 4 - 
25898331375g 2 g^g 4 + 47959232650g^g^ - 33861927750Q^Q|Q|+ 
7454763225g^g^ - 184988475g 2 , g 4 + 11891917800Q^Q|Q 5 - 

36io54638oog^g^g 4 g 5 + 27953667ooog^g^g 5 - 44577i62oog^g 3 g^g 5 + 

6285855240g^g^g^ - 6471756144g^g2g 4 g2 + 565259688g^g^g^+ 
380939328g^g 3 g^ - 4375251160g^g 6 + 10317018600g^g 4 g 6 - 

5ii38i332og 2 , g 3 g 4 g 6 + 23888844og 2 1 g 4 g 6 - 320387i552g^g^g 5 g 6 + 

1758685824g^g 3 g 4 g 5 g 6 - 88566912Q£Q|!Q 6 + 335466432g 2 1 g 3 g6- 
55073088g^g 4 g^ + 1351294560g 2 t g 3 g 7 - 2341442880g^g^g 4 g 7 + 

62654256og^g 3 g^g 7 + 6i9520832g 2 1 g 3 ; g 5 g 7 - i23524352g^g 4 g 5 g 7 - 

96574464g^g 3 g 6 g 7 + 4048704g^ - 341160120g^g^g 8 + 

3862ioi6og^g 4 g 8 - 3083724og 2 / g 4 g 8 - 78073632g^g 3 g 5 g 8 + 
58484i6g^g 6 g 8 + 704i552og 2 1 g 3 g 9 - 4342464og^g 3 g 4 g 9 + 52557i2g®g 5 g 9 )- 

(n + l/2) 6 i 

1 -72576Qi 2 Ql + 1886976QnQ|g 3 - 



300987187200g^ 2 v^g2" 
22135680g 10 g2g 3 + 27463538375g^° + 2903040gi g^g 4 - 
141448688325g 2 g^g 4 + 240655765350g 2 g6g2 _ 152907158250g^g^+ 
28724479875g^g^g^ - 413669025g 2 , g 4 + 59058073200g^g^g 5 - 

i642642092oog^g^g 4 g 5 + ii36546964oog 2 t g 3 g 4 g 5 - i5i663428oog 2 , g 3 g 4 g 5 + 
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26061194880Q4Q4Q2 _ 23876233920Q^q 4 q2 + i767i89312Q^g2 + 

1292433408Q^Q 3 <55 - 18902165520Q^Q 6 + 40256773200Q^Q^Q 4 Q 6 - 
1711697A000QlQlQlQ 6 + 483582960Q^Q^g 6 - 11384150400Q^Q 5 g 6 + 
5285056896Q^Q 3 Q4<25<26 - 246903552Q^Q 2 ,Q 6 + 992779200Q^Q^g2_ 
101860416Q^Q 4 Q 2 6 + 4966859520Q^Q^Q 7 - 766i606400Q^Q 3 i Q 4 Q7+ 
1683037440g^g 3 <5lQ7 + 186157 U00Q 6 2 Q 2 3 Q 5 Q 7 - 316U1056Q 7 2 Q 4 Q 5 Q 7 - 
235U62A0Q 7 2 Q 3 Q 6 Q 7 + 88957UQ 8 2 Q 2 7 - 1042372800Q^Q 8 + 
1016789760Q^Q^Q 4 Qs - 52436 160Q^Q 8 - 189060480Q^Q 3 Q 5 Q 8 + 
9217152Q*Q 6 Q 8 + 175190400Q^Q^Q 9 - 87816960Q^Q 3 Q 4 Q 9 + 10378368Q^Q 5 Q 9 ) (16) 
All six WKB corrections printed in MATEMATICA are available from the author in 



6 Appendix 2: QNMs of a 4- dimensional Schwarzshild 
black hole 

"The potential" Q(x) in case of a Schwarzshild black hole has the form 



where s = corresponds to scalar perturbations, s = 1/2 - neutrino perturbations, s = 1 
- electromagnetic perturbations, s = 3/2 - gravitino perturbations, s = 2 - gravitational 
perturbations. The quasinormal frequencies at 3th and 6th WKB orders and in compar- 
ison with numerical results [HS! are presented in the table I. 
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(17) 
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s - 


= 




numerical 


3th order WKB 


6th order WKB 


I = 0, 


n = 





0.1105 


- 0.1049i 


0.1046 - 


0.1152vj 


0.1105 


- 0.1008; 


I = 1, 


n = 





n 9Q9Q 


n no 77.?' 
— u.uy ( 1 1 


n 9Qi i 


u.uyoui 


n 9Q9Q 

u.zyzy 


n no 78i 
— u.uy ( oi 


/ — 1 
t — 1, 


Tl = 


i 


0.2645 


- 0.3063i 


0.2622 - 


0.3074; 


0.2645 


- 0.3065; 




71 — 


o 


0.4836 


- 0.0968i 


0.4832 - 


0.0968; 


0.4836 


- 0.0968; 


/ — 9 




1 


0.4639 


- 0.2956i 


0.4632 - 


0.2958i 


0.4638 


- 0.2956; 


/ — 9 




9 


0.4305 


- 0.5086i 


0.4317- 


0.5034; 


0.4304 


- 0.5087; 


s = 


1 /O 

1/2 




numerical 


3th order WKB 


6th order WKB 


I = 1, 


n = 









0.2803 - 


0.0969; 


0.2822 


- 0.0967; 


7 1 

' — 1) 


n = 


1 






0.2500 - 


0.3049; 


0.2525 


- 0.3040; 


7 O 

/ = 2, 


n = 









0.4768 - 


0.9639; 


0.4772 


- 0.0963; 


7 O 

/ = 2, 


n = 


1 






0.4565 - 


0.2947; 


0.4571 


- 0.2945; 


7 O 

/ = 2, 


n = 


2 






0.4244 - 


0.5016; 


0.4231 


- 0.5070; 


1 — 3 

l — O, 


TL — 


n 

U 






0.6706 - 


0.0963; 


0.6708 


- 0.0963; 


/ = 3, 


Tl — 


1 






0.6557- 


0.2917; 


0.6560 


- 0.2917; 


/ — 3 


in — 
lb — 


o 






0.6299 - 


0.4931; 


0.6286 


- 0.4950; 


/ — 3 




o 






0.5970 - 


0.6997; 


0.5932 


- 0.7102; 


5 = 


= 1 




numerical 


3th order WKB 


6th order WKB 


7 1 

/ = 1, 


n = 





0.2483 


- 0.0925i 


0.2459 - 


0.0931; 


0.2482 


- 0.0926; 


7 1 

/ = 1, 


n = 


1 


0.2145 


- 0.2937i 


0.2113 - 


0.2958; 


0.2143 


- 0.2941; 


7 C\ 

I = 2, 


n = 





0.4576 


- 0.0950i 


0.4571 - 


0.0951; 


0.4576 


- 0.0950; 


7 C\ 

I = 2, 


n = 


i 

1 


0.4365 


- 0.2907; 


0.4358 - 


0.2910; 


0.4365 


- 0.2907; 


1 o 

' — 2, 


n = 


2 


0.4012 


- 0.5016i 


0.4023 - 


0.4959; 


0.4009 


-0.5017; 


/ — 3 


Tl 


n 

u 


0.6569 


- 0.0956i 


0.6567- 


0.0956; 


0.6569 


- 0.0956; 


/ = 3, 


■ = 


1 


0.6417 


- 0.2897i 


0.6415 - 


0.2898; 


0.6417 


- 0.2897; 


/ — 3 


Tl 


9 


0.6138 


- 0.4921i 


0.6151 - 


0.4901; 


0.6138 


- 0.4921; 


/ — 3 


11 


Q 
O 


0.5779 


- 0.7063i 


0.5814- 


0.6955; 


0.5775 


- 0.7065; 


S = 


3/2 




numerical 


3th order WKB 


6th order WKB 


7 1 

' — 1) 


Tl = 









0.1817- 


0.0866; 


0.1739 - 


- 0.08357; 


7 1 

' — 1) 


Tl = 


1 






0.1354- 


0.2812; 


0.1198 


- 0.2813; 


7 O 

/ = 2, 


Tl = 









0.4231 - 


- 0.926; 


0.4236 


- 0.0925; 


7 O 

/ = 2, 


Tl = 


1 






0.4000 - 


0.2842; 


0.4007 


- 0.2838; 


/ = 2, 


Tl = 


2 






0.3636 - 


0.4853; 


0.3618 


- 0.4919; 


7 — 3 

l — O, 


Tl = 


n 

U 






0.6332 - 


0.0945; 


0.6333 


- 0.0944; 


/ = 3, 


Tl 


1 






0.6173- 


0.2864; 


0.6175 


- 0.2863; 


/ — 3 


/ v — 


9 






0.5898 - 


0.4846; 


0.5884 


- 0.4868; 


/ — 3 


/ v — 








0.5547 - 


0.6882; 


0.5505 


- 0.7000; 


s = 


= 2 




numerical 


3th order WKB 


6th order WKB 


7 O 

/ = 2, 


Tl = 





0.3737 


- 0.0890i 


0.3732 - 


0.0892; 


0.3736 


- 0.0890; 


/ = 2, 


Tl = 


1 


0.3467 


- 0.2739i 


0.3460 - 


0.2749; 


0.3463 


- 0.2735; 


/ = 2, 


Tl = 


2 


0.3011 


- 0.4783i 


0.3029 - 


0.4711; 


0.2985 


- 0.4776; 


I = 3, 


Tl = 





0.5994 


- 0.0927i 


0.5993 - 


0.0927; 


0.5994 


- 0.0927; 


7 O 

/ = 3, 


Tl = 


1 


0.5826 


-0.2813i 


0.5824 - 


0.2814; 


0.5826 


- 0.2813; 


/ = 3, 


Tl = 


2 


0.5517 


- 0.479H 


0.5532 - 


0.4767; 


0.5516 


- 0.4790; 


7 O 

I = 3, 


Tl = 


3 


0.5120 


- 0.6903i 


0.5157- 


0.6774; 


0.5111 


- 0.6905; 


7 A 

I = 4, 


Tl = 





0.8092 


- 0.0942i 


0.8091 - 


0.0942; 


0.8092 


- 0.0942; 


7 — 4 


Tl = 


i 
i 


0.7966 


- 0.2843i 


0.7965 - 


0.2844; 


0.7966 


- 0.2843; 


I = 4, 


Tl = 


2 


0.7727 


- 0.4799i 


a.?736 - 


0.4790; 


0.7727 


- 0.4799; 


Z = 4, 


Tl = 


3 


0.7398 


- 0.6839i 


0.7433 - 


0.6783; 


0.7397 


- 0.6839; 


Z = 4, 


Tl = 


4 


0.7015 


- 0.8982i 


0.7072 - 


0.8813; 


0.7006 


- 0.8985; 
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